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Abstract 

We study the arrival time distribution of overdamped particles driven by a constant force in a 
piecewise linear random potential which generates the dichotomous random force. Our approach 
is based on the path integral representation of the probability density of the arrival time. We 
explicitly calculate the path integral for a special case of dichotomous disorder and use the cor- 
responding characteristic function to derive prominent properties of the arrival time probability 
density. Specifically, we establish the scaling properties of the central moments, analyze the be- 
havior of the probability density for short, long, and intermediate distances. In order to quantify 
the deviation of the arrival time distribution from a Gaussian shape, we evaluate the skewness and 
the kurtosis. 

PACS numbers: 05.40.-a, 05.10.Gg, 05.60.-k 
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I. INTRODUCTION 



The overdamped equation of motion of a classical particle in one dimension presents a 
simple yet very useful model for the study of many physical, biological, economical, and 
other systems. Depending on the character of the force field acting on the particle, this 
equation provides a basis for describing different phenomena in these systems. Specifically, 
if the force field contains the noise terms arising from the influence of the environment, 
it describes a large variety of noise-induced phenomena including noise-induced transitions 
g directed transport B, and stochastic _ g to name on! y a few. It snonld be 
noted that in some cases, especially within the white-noise approximation, the statistical 
properties of the solution of this equation can be obtained analytically. 

If the environment is disordered the force field contains also the random functions of the 
spatial variable. In this case the overdamped dynamics represents both the noise-induced 
and disorder-induced effects and it can exhibit as well anomalous behavior even in the 
simplest situation of additive white noise [4j|. For the latter situation a number of exact 
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results was obtained for Sinai disorder 
and also for some special cases of non-Gaussian disorder 

When the noise terms produced by a stochastic environment become negligible, the over- 
damped equation of motion accounts solely for effects of quenched disorder. This equation 
effectively describes, e.g., the transport of particles in deterministic ratchets with quenched 
disorder Ji], 
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2 if ] and can be used for the study of the dynamics of localized structures 



like domain walls in random magnets and vortices in type-II superconductors. Although 
temporal noise terms are absent, there are only very few exact results available. Therefore, 
in order to fill this gap, we have examined the following dimensionless equation of motion 
for an overdamped particle 221 ]: 

X t = f + g(X t ). (1.1) 

Here, X t denotes the particle coordinate that satisfies the initial condition Xq = 0, / (> 0) 
is a constant force, and g(x) = —dll(x)/dx = ±g is a dichotomous random force generated 
by a piecewise linear random potential U(x), see Fig. 1. It is assumed that the random 
intervals Sj of a linearly varying U(x) are statistically independent and distributed with the 
same (exponential) probability density p(s). Moreover, we assume that the conditions / > g 
and g(+0) = —g are imposed. 
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Equation (11. ip is of minimal form that accounts for the effects of quenched disorder 
on the overdamped motion of driven particles. Its main advantage is that many of the 
statistical properties of X t can be described analytically in full detail. Nevertheless, if the 
odd and even intervals Sj are distributed with different exponential densities (in this case 
the exact results exist as well), then Eq. (1.1) can be used also for studying a number 
of important physical issues. Specifically, this equation constitutes a basis for describing 
the adiabatic transport of particles in randomly perturbed one-dimensional channels and 
presents a simple model for studying the low-temperature dynamics of charge carriers and 
localized structures in randomly layered media. In addition to the listed examples, we point 
out at rather unexpected application of Eq. (1.1) in astrophysics. Namely, if the clouds in 
interstellar space are distributed uniformly then the distances between them are distributed 
with an exponential distribution. In this case, assuming that the light velocity in the clouds 
is the same and the sizes of clouds are distributed with an exponential distribution, Eq. (1.1) 
can be used for studying the statistical properties of distances that pass the light emitted 
by a star in different directions. 



In [22|] we derived the probability density of the solution of Eq. ( 11.11) and investigated 
explicitly its time evolution. In contrast, in this work we focus on the statistical properties 
of the arrival time for the particles governed by Eq. (11.11) . The paper is structured as follows. 
In Sec. II, we derive the path integral representation for the probability density of the arrival 
time. The characteristic function of the arrival time is determined in Sec. III. In Sec. IV, we 
calculate the moments of the arrival time and study their asymptotic and scaling behavior. 
The basic properties of the arrival time probability density are studied in Sec. V, both 
analytically and numerically. We summarize our findings in Sec VI. Some technical details 
of our calculations are deferred to the Appendix. 

II. PATH INTEGRAL REPRESENTATION OF THE ARRIVAL TIME PROBA- 
BILITY DENSITY 

According to Eq. (11. 1ft . the arrival time t x , i.e., the time that a particle spends moving 
from the origin to a position x(> 0), is given by the integral expression 

f x dx 
Jo f + g{x) 
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This time depends on the random function g(x) and thus presents a random quantity. The 
probability density P x (t) that t x = t for a fixed coordinate x, i.e., the probability density of 
the arrival time, is defined in the well-known way as 

P x (t) = (5(t-t x )), (2.2) 

where the angular brackets denote an averaging over the sample paths of g(x), and S(t — t x ) 
is the Dirac 5 function. 

To obtain the explicit form of P x (t) we use a path integral approach. Because of the 
dichotomous character of the random function g(x), it is convenient to present the probability 
density in terms of the partial densities 

oo 

ft(«) = E p i B) (*), (2-3) 

n=0 

where p( n '(t) is the probability density that for the sample paths of g(x) which are under- 
going n changes of the sign on the interval (0, x) the condition t x = t holds. For a given 
n (> 1) the solution of Eq. (II. ip can be written in the form 

4 n) = E S / + Vi ( 2 - 4 ) 

with s n+ i e (0, s n+ i). On the other hand, because g(x) = (— l) 3 g if x belongs to the interval 
Sj, flUD yields 

An) _ y- s j , Co pil 

" "fe/+(-l)^ + / + (-l)^' ^ j 

Setting = x and replacing s n+ i by x — Sj=i Sj, the result (12.51) can be recast to 

^ = T^-f^¥W- ir + ( - 1)% - (2 ' 6) 

Let us next introduce the probability p(sj)dsj that the j-th jump of g(x) occurs in 
the interval dsj and also the probability ff° p(s)ds that the distance between the nearest- 
neighbor jumps exceeds I. Then, the probability dW n (x) that the function g(x) on the 
interval (0, x) experiences n jumps in the intervals dsj (j = 1, . . . , n) assumes the form 

/■oo n 

dW n {x)= p{s)dsl[p(s j )ds j . (2.7) 



3=1 
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Because s n+ i > 0, the positive variables of integration, Sj, must satisfy the condition 
J21=i s j < x - Denoting by £l n (x) a region in the n-dimensional space of these variables, 
being defined by the aforementioned condition, we obtain 

PW(t)= / r 5(t-tP)dW n (x). (2.8) 

Finally, taking into account that = x/(f — g) is the arrival time at n = and Wo 0*0 = 
J£° p(s)ds is the total probability of those sample paths of g(x) which do not change sign 
on the interval (0,x), we end up with the following path integral representation for the 
probability density of the arrival time: 

oo „ 

P x (t) = 6{t - tf)W {x) + £ / 5(t- t^)dW n (x). (2.9) 

This form of the arrival time probability density is rather general, but possesses a rather 
complex mathematical structure. Based on (12. 9p . however, we arrive at two conclusions that 
are valid for an arbitrary probability density p(s): (i) P x (t) at a fixed x is concentrated on 
the interval [x/(f + g),x/(f — g)], and (ii) P x (t) is properly normalized, i.e., / °° P x (t)dt = 1. 
Indeed, since minij^ = x/(f + g) and max^"' — x/(f — g), we have 5(t — t^ 1 ') = and 
so P x {t) = if t ^ [x/(f + g),x/(f — g)\. To prove the second assertion, we first note 
that, according to (ED, / °° = W (a;) + E£=i W„(x), where W n (x) = dW n (:r) 

is the probability that the function g(x) has undergone n jumps in the interval (0, x). Next, 
introducing the quantities S n (x) = TT?-i p(sj)dsj, we find the representations Wo(x) = 
1 — S±(x) and W n (:r) = S^x) — S n+ i(x), see also [22J. Finally, taking into account that 
Soo(x) = and Y^Li W n (x) = Si(x), we assure that the normalization condition holds true 
for an arbitrary p(s). 

III. CHARACTERISTIC FUNCTION OF THE ARRIVAL TIME 

By use of the integral formula for the 5 function, 

5(t - 4 n) ) = r e-^'^du, (3.1) 

27T J-oo 

we can rewrite the probability density (I2.9P in the form of a Fourier integral, i.e., 



Px(t) = ^-I M^)e- tuJt dw. (3.2) 

Z7T J-oo 
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According to Eqs. (I2.9p - fl3.2l) . the characteristic function <p x (uj), which determines all the 
statistical properties of the arrival time t x , is obtained as 

<f) x ^)=e luJti0) Wo(x)+Muj), (3.3) 

where 

Uu) = £ / e^" dW n (x). (3.4) 

In the general case of an arbitrary p(s), the characteristic function has a complex structure 
involving an integration over the n-dimensional domain Q n (x) and a summation over all n. 
Remarkably, however, 4> x (uj) can be expressed in terms of elementary functions if the random 
intervals Sj are exponentially distributed, i.e., if p(s) = \e~ Xs , where A -1 is the average length 
of Sj. In this case the probability (12.71) becomes 

n 

dW n (x) = e~ Xx X n J] d Sj , (3.5) 
3=1 



and (13.41) reduces to 



00 r m n 

-Xx \ y \n I Aujt 



_[fj) = e J2 Xn , lidsj. (3.6) 

n=l J U n (x) - =1 

For the calculation of (j) x (u)) it is convenient to transform the right-hand side of Eq. (13. 6p 
into a form with separate integrations over the variables Sj. To this end, we use an approach 



22j based on the integral representation of the step function 



1 r+00 e (i*+ri)y 1 if y > 

■dK={ (3.7) 

if y < 



2tt J -00 in + rj 



which is valid for arbitrary 77 > 0. Applying (13. 7p to (13.61) and setting y = x — J2]=i s j, one 
obtains the desired result: 



^ W ) = ^VA"/ dK- / ^ 

2n ^ J -00 ik+ r\ Jo Jo 



xe -(^)E;=i^- f[ dSj . (3. 8 ) 

Then, using the identity J2%Li a n = Em=i[ ffl 2m-i + a 2m] and taking into account that, 
according to (12. 6p . 

y.(2m-l) _ X , 2 # c 



we reduce the formula (13. 8p to the form 



g-Arr .qo ^uox _co 



— E / m ~Vo)-H^i)e W(/+9) 

-co Vft 1 '- 



Here, 



{Rev k > 0, k = 0,1,2) and 



2vr 7-oo u m=1 
+ J m (i/ )/ m (^)e W(/ ~ 9) ]^. 

*K) = r p (s)e-^- x >ds = - 
Jo v k 



= m + ?7 + -d fc 



(3.10) 



(3.11) 



(3.12) 



P-9 2 

with 5o = 0, Si — — 1, and S 2 — 1. We note that the right-hand side of (13.101) contains an 
arbitrary positive parameter 77. According to the definition (13.41) . however, the left-hand side 
of (13.101) does not depend on rj. This implies that the final result of evaluating the series 
and integral in (13.101) does not depend on rj as well. Therefore, for auxiliary manipulations 
we may choose a most convenient value for this parameter. 

From this point of view, it is reasonable to choose rj > A. This is so because in this case 
1-^(^)1 < 1 an d the series in (13.101) can be easily evaluated: 

L - -{Uo)l -{Vl) = — 

m=l 



£ J m - 1 (z/ )/ m (^i) 

1=1 

00 

£ r>o)r> 2 ) 



m=l 



Substituting (13.131) into (13.101) and using that W (x) = e and 



VqV x - A 2 ' 

A 2 

v Q v 2 - A 2 ' 

—Ax 



(3.13) 



(0) 



e 1 ^ Wo (x) 



2tt 



00 Z/q 



(3.14) 



from ( 13.31) we obtain 



-Ax 



2tt 



dn. 



(3.15) 



z/ z/i - A 2 z/ ^2 - A 2 
Upon calculating the integrals in (13.151) (the details are given in the Appendix) we find 
a remarkably simple expression for the characteristic function of the arrival time 



'Xx(l-ivffg) 



cosh 



— V 



+ 



1 +iv 



sinh 



(3.16) 



where 



' KP-9 2 Y (3,17) 

We note that (f> x {w), being the characteristic function, satisfies the conditions |0 a; (a')| < 1, 
4> x (0) = 1 and 4>x(— oj) = <j>%{w), where the asterisk indicates the complex conjugation. 
Equation (13.161) is our main result which allows us to study the essential properties of the 
arrival time analytically. 



IV. MOMENTS OF THE ARRIVAL TIME 



The moments of the arrival time are defined in the usual way as (t™) = /f^ t m P x (t)dt, 
and can be deduced through the characteristic function as follows: 

1 d m 



uj=0 



i m du m 

According to (13.161) and (|4.1|) . the first moment, i.e., the mean arrival time, emerges as 

1 



(4.1) 



2Xfx + g - ge 



-2 Ax 



(4.2) 



At small distances from the origin, when Xx 1, the formula (14. 2p yields (t x ) = x/(f — g). 
This result is expected: The total probability of those sample paths of g(x) which do not 
change the sign on the interval (0, x) tends to 1 as Xx — > 0, and thus the average particle 
velocity tends to f—g. In the other limiting case, when Xx ^> 1, the formula (14.2ft approaches 
(tx) — f x /{f 2 ~ 9 2 )- This result is corroborated by the fact that the long-time asymptotic 
of the average particle velocity equals (f 2 — g 2 )/ f [22] . 

The moments of higher order can also be calculated straightforwardly. In particular, for 
the second moment we obtain 

1 



(tl) 



2X 2 (f 2 -g 2 ) 2 
- g (2Xfx-g)e- 2Xx 



2X 2 f 2 x 2 + 2Xg(f + g)x-g< 



(4.3) 



The central moments, ((t x — (t x )) m ), can be determined from the finite series 



((*.-<or>=E(- i ) tB " i csi.^><*»> 

j=0 



m—l 



(4.4) 
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where C 3 m is the binomial coefficient, or, alternatively, by the formula 

((t x - (t x )) m ) = ; — <j)M e _MU ■ 

u x x " 1 i m du) m ^=o 

Specifically, using either of these definitions, the variance of the arrival time, a 2 = ((t 

(tx)) 2 ), can be written in the form 



(4.5) 



9 Z 



AXx - 3 + 4e~ 



At short distances, when Xx <C 1, the variance reduces to 



2Xx _ g-4Ax 



KP-9 2 ) 2 

and at long distances, when Xx 3> 1, it reduces to 



x 



a' 



if 



-x. 



(4.6) 



(4.7) 



(4.8) 



x KP-9 2 ) 2 ' 

Moreover, the central moments of the arrival time possess interesting scaling properties. 
Namely, using (14.51) with (13.161) and (I4.2p . one obtains 

9 



where 



^m(Ax) 



KP-9 2 ) 



^m(Xx), 



Jm 

- Xr ' <!>(.;. A.r) 



2=0 



(4.9) 



(4.10) 



is a function of the single variable Ax, and 

/ 1 — e~ 2Xx 
&(z, Xx) = exp ( — z 



cosh ( AxVl + z 2 



+ 



l + z 

ViTz 



: sinh 



(4.11) 



Thus, the central moments exhibit a universal dependence on /, g and A, i.e., ((t x — (t x )) m ) oc 
[g/X(p - g 2 )] m . 



V. PROPERTIES OF THE ARRIVAL TIME PROBABILITY DENSITY 

As it follows from (I3.16p . the characteristic function tends to e~ Xx + wx l^-9) as \u\ — > oo. 
According to (I3.2p . this suggests that the probability density of the arrival time contains the 
(^-singular contribution, i.e., 

Pxit) =s(t- j^A e- Xx + P x (t), (5.1) 



where 

Px(t) = — / <t> x (u)e du (5.2) 

Z7T J-oo 

denotes the regular part of P x {t) and 

4> x {u) = <j> x {u) - e -^+W(/-9). (5.3) 

Because the intensity of the 5-singular part decreases exponentially with increasing x, its 
contribution to P x (t) plays a crucial role only at short distances from the origin. The regular 
part rules the behavior of P x (t) at longer distances. 



A. Behavior at short distances 



At Ax 1 we can obtain the probability density P x (t) in a simple way, without the need 
to evaluate the integral in (15. 2p . To this end, we first note that the 5-singular part of P x (t) 
is formed by those sample paths of g(x) that do not change the sign on the interval (0,x). 
Accordingly, only the sample paths which have at least one change of the sign on this interval 
do contribute to the regular part P x (t). For small values of Ax repeated changes of the sign 
are unlikely. Therefore, in order to determine P x (t), we consider the sample paths with a 
single change of the sign. In this case the probability P x (t)dt is equal to dW\(x) = \e~ Xx ds\ 
and, because t = t^' = x/(f + g) + 2gsi/ (f 2 — g 2 ), the relation P x (t)dt = dWi(x) at Ax 1 
then yields 

P x (t) = (5.4) 

Finally, substituting (15 .4p into (15.11) . we find the probability density of the arrival time at 
Ax < 1: 

P x (t) =s(t- j^j (1 - Ax) + HP 2g g2) . (5.5) 

At first sight, this result may come as a surprise because the regular part of the probability 
density does not depend explicitly on x and t. It should be stressed, however, that the 
formula (15.41) is derived under the condition that t G [t m in,^max], where t m i n = x / '(/ + g) and 
^max = x /(f — g) [ we recall that P x (t) = if t ^ [t m in, ^max]]- This means that P x (t) depends 
on x and t implicitly leading to the broadening of P x (t) if x increases. At the starting point 
x = we have t min = t max = 0, therefore Po(t) = and, in accordance with the condition 
to = 0, Po(t) = S(t). We note also that the normalization condition //™ P x (t)dt = 1, which 
holds true also for (15.51) . further corroborates the validity of (15 .4p . 

10 



B. Behavior at long distances 



To study the long-distance behavior of the probability density P x {t), it is convenient to 
introduce the new time variable r = (t — (t x ))/cr x . The corresponding scaled probability 
density V x {t) is expressed through P x (t) as V x {j) = cr x P x ({t x ) + a x r) and, according to 
(13. 2p . it can be written in the form 

V x {r) = — / M»M e-^^-'^dfi. (5.6) 

Using ( 13. 16[) . ( 14.21) and ( 14. 61) . the characteristic function of V x (t) then reads 

<f>M(T x ) e-^l°* = e- Al '$(z/i^ 2 1/2 (Ax), Ax), (5.7) 

where the function $(z, Ax) is defined by Eq. (14. lip and 

^ 2 (Ax) = Ax - ~ + e~ 2Xx - ~ e~ 4Xx . (5.8) 

Because the characteristic function (15. 7p depends only on Ax and the integration variable 
/i, the scaled probability density ( 15. 6p possesses the remarkable property that V x [j) is a 
function of Ax and r which depends neither on the external force / nor on the amplitude g 
of the dichotomous random force g(x). 

In the case of long distances, if Ax 3> 1, the characteristic function (15. 7p at /x 4 Ax can 
be approximated by the two terms of its expansion: 

faitiM e"^>/- = e^ 2 / 2 ( 1 - (5.9) 

\ 8Ax J 

Substituting (15. 9p into (15.71) and calculating the integrals, we find the two terms of expression 
of the scaled probability density 

which is valid if Ax ^> max(l,r 4 ). Thus, in accordance with the central limit theorem 
of probability theory (see, e.g., Ref. 231]). the limiting probability density approaches a 
Gaussian form, i.e., 'P 00 (t) = (27r)~ 1 / 2 e~' r2//2 , and V x {r) — Voo(r) oc (Ax) -1 as Ax — ► oo. 



C. Numerical verification 



Our numerical calculations pursue two goals, namely (i) to verify the analytical findings 
and (ii) to illustrate and visualize the obtained findings. The former is achieved by compar- 
ison of the probability density (15.11) with that derived from the numerical simulation of the 
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arrival time (12.11) . We use the Maple package for calculating the Fourier integral in (I5.2p 
and employ the histogram procedure to numerically evaluate the probability density. In 
short, this procedure consists of successive generations of random intervals Sj according to 
the exponential distribution and evaluating the arrival time to a fixed position x for different 
realizations of random intervals. The probability density is then presented as the histogram 
of arrival times of the particle. For further details about this procedure we refer the inter- 
ested reader to Ref. [22] where a similar approach was used for the numerical evaluation of 
the probability density of the particle position at a fixed time t. In doing so, we made sure 
that the simulated probability density function is in perfect agreement with the theoretical 
one, see Fig. 2. 

Figure 3 illustrates the short- distance behavior of the probability density fl5.ll) . As can 
be seen in Fig. 3(a), at very small values of x the probability density of the arrival time is 
described by the approximate formula (15.51) . In accordance with the assumption made in its 
derivation, this suggests that the sample paths of g(x) which on the interval (0, x) have more 
than one change of the sign are responsible for the explicit dependence of P x (t) on t. If x is 
not too small, i.e., the total probability of these sample paths is small but non-zero, then, as 
shown in Fig. 3(b), P x (t) is an almost linear function of t. With increasing of x the role of 
these sample paths becomes increasingly important: The function P x (t) becomes nonlinear, 
assumes a unimodal form, and eventually approaches a Gaussian shape, see Fig. 4. 



D. Skewness and kurtosis 



In order to quantitatively describe the difference between the arrival time probability 
density and a Gaussian density with identical mean and variance as P x (t), we calculate the 
skewness 

, (J) = (g- - <«.»'> (5.11) 

x 

that characterizes the degree of asymmetry of P x (t), and as well the kurtosis 

k(x) = ^ tx ~ f * ))4) - 3 (5.12) 

that characterizes the degree of peakedness of P x (t). Because s(x) = and k(x) = if the 
arrival time t x follows a Gaussian distribution, one can consider the skewness and kurtosis 
as appropriate measures of deviation of the arrival time distribution from a Gaussian shape. 
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Using the representation (14. 9 j) for the central moments, from the definitions (15.111) and (I5.12p 
we obtain 

six) - * 3{XX) kix) - ^ 4(AX) - 3 (5 13) 

i.e., s(x) and k(x) are universal functions of the single variable Ax, see Fig. 5. Calculating 
^3 (Ax) and \l/ 4 (Ax) and taking into account (I5.8p . we find an explicit expression for the 
skewness, 

= ~(4Ax-3 + 4e- 2A *-e- 4A *) 3 / 2 ^ + (3 

- 12Xx)e- 2Xx - 6e~ iXx + e~ 6Xx ], (5.14) 

and for the kurtosis, 

6 

kix) = — — —. .. Nn 13 — 8Ax 

V ; (4Ax - 3 + Ae~ 2Xx - e ~ 4Xx ) 2 1 

- 8(1 + AXx)e- 2Xx - 4(3 - 4Xx)e~ 4Xx 

+ 8e- 6Xx - e- 8Xx }. (5.15) 

The formulas (I5.14j) and (15.151) yield in leading order of Ax the following relations: s(x) = 
-(3 v / 3/4)(Ax)- 1 / 2 and k(x) = (9/5)(Ax)~ 1 at Ax < 1, and s(x) = -(l/2)(Ax)- 3 / 2 and 
k{x) = — 3(Ax) _1 at Ax 3> 1. These results clearly evidence that the arrival time probability 
density P x (t) distinctly differs from a Gaussian density at short distances and approaches 
this Gaussian shape at long distances. Moreover, since \k(x)/s(x)\ — > oo as both, Ax — > 
and Ax — >• oo, the kurtosis can be considered as a unique measure of non-Gaussianity of 
P x (t). We note that, because of the condition s(x) < 0, the left tail of P x (t) is always 
heavier than the right tail. Also, P x (t) is more peaked compared to the Gaussian density at 
distances where k(x) > 0, and is more flattened at distances where k(x) < 0. 



VI. CONCLUSIONS 



We applied the path integral approach to calculate the characteristic function of the 
arrival time for overdamped particles driven by a constant bias in a piecewise linear ran- 
dom potential producing a dichotomous random force with exponentially distributed spatial 
intervals. Using the characteristic function, we derived the moments of the arrival time, 
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established universal scaling properties of the central moments, and demonstrated that the 
arrival time probability density P x (t) contains both a 5-singular contribution and a regular 
part. While the (^-singular part, whose weight decreases exponentially with increasing x, 
plays the main role at short distances, the regular part of P x (t) dominates at large distances 
x. 

At very small distances the regular part is defined by the sample paths with only one 
change of the sign on the interval (0, x) and in this case its value does not depend on x. Upon 
increasing x the contribution of other sample paths leads to the transformation of this part 
of P x (t) into an almost linear function oft, and subsequently into unimodal form, and finally, 
at x — > oo, it tends to a Gaussian density as x^ 1 . Moreover, in order to characterize the 
difference of the arrival time probability density from the Gaussian density, we calculated the 
skewness and kurtosis. The function P x (t) is more peaked in comparison with the Gaussian 
density at small x and is more flattened at large x. 
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APPENDIX: DERIVATION OF THE CHARACTERISTIC FUNCTION 



For calculating the integrals in (13. 15j) . 



r = J-f 

27T J-oo 
1 



\ e v x 
\2 1 



2vr 



we use the method of contour integration 



u 2 e 



U X 



u u 2 - A 



2 dtij 



(A.l) 



24| . According to (I3.12p and (13.171) , the integrands 



in (IA.ll) . R(k) = Xe u ° x / (vqV\ — A 2 ) and S(k) = u 2 e v ° x / {v§v 2 — A 2 ), can be written in the form 

A 



R(k) 



14 



S(k) 



in + i] + HXv 



(A.2) 



where 



K i,2 = irj ± iXyi — v 2 + Xu, 
K 3,4 — if] i iXy/l — v 2 — Xu. 



(A.3) 



The formulas ( 1A.2I) exhibit that both i?(ft) and £(«;) as functions of the complex variable 
k have two poles of the first order at k = K\p and k = ^3,4, respectively. If 77 > A then all 
poles are located in the upper half plane of the complex k plane, and the residue theorem 
yields 



Y = z[Resityci) +Resi?(« 2 )], 
Z = i[Res S(ks) + Res S(k,4)]. 



(A.4) 



Since the residues in (I A. 4 1) are defined as Res-R^i^) = lim K _ >Kl2 (K — ki i2 )R(k) and 
Res5 , (K 3i4 ) = lim K ^ K34 (K 



K3 i 4)S , (k), from ( 1A.2I) and ( 1A.3I) we obtain 



F = , sinh (XxV 1 — i/ 2 ) e' 1 



sinh 



+ cosh 

Finally, substituting ( 1A.5I) into the formula 



-iXxu 



(A.5) 



u; = e 



-A:r 



y e «W(/+9) _|_ Z e iux/(f-9) 



(A.6) 



which follows from (I3.15P and flA.ll) . we get the desired characteristic function (13.161) . 
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FIG. 1: Schematic representation of (a) the piecewise linear random potential U(x) and (b) the 
corresponding dichotomous random force g(x) = —dU(x)/dx as functions of the coordinate x. 
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FIG. 2: (Color online) Theoretical and simulated probability density of the arrival time for x = 2.4. 
The solid line (red online) and histogram represent the analytical result (15. ip and the numerical 
simulation of the arrival time (|2.ip . respectively. The parameters of the force field are chosen as 
/ = 1, g = 0.3, and A = 1. The vertical arrow depicts the ^-singular contribution to P x (t). 
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FIG. 3: Short-distance behavior of the probability density of the arrival time P x (t) at (a) x = 10 -3 
and (b) x = 0.2. The other parameters are the same as those in Fig. 2. 




FIG. 4: (Color online) Plots of the probability density of the arrival time for different distances 
from the origin of the coordinate system. The vertical surfaces (green online) depict the regular 
part of P x (t). In order to visually demonstrate that the intensity of the (^-singular part of P x (t) 
exponentially decreases with x, we depicted the length of the vertical arrows (red online) in the 
form 0.8 + 2e~ Xx . For the convenience of comparison, the force field characteristics are chosen as 
in Figs. 2 and 3. 
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FIG. 5: Plots of the skewness s(x) (dashed line) and the kurtosis k(x) (solid line) of the arrival 
time probability density vs the normalized distance Ax. 
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